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Recent thermal conductivity data on the heavy fermion superconductor UPtg have been interpreted 
as offering support for an E2u model of the order parameter as opposed to an Eig model. In this 
paper, we analyze this issue from a theoretical standpoint including the detailed effects of Fermi 
surface and gap anisotropy. Our conclusion is that although current data put strong constraints on 
the gap anisotropy, they cannot definitively distinguish between these two models. Measurements 
on samples of varying quality could be decisive in this regard, however. 



PACS Numbers: 74.70Tx,74.25Fy 

Well over a decade after the discovery of heavy fermion 
superconductivity, the pairing mechanism and even the 
order parameter symmetry in these compounds remain 
controversial. Early suggestions of pairing in an uncon- 
ventional superconducting state, based primarily on anal- 
ysis of transverse ultrasound measurements in UPt^, 0] 
were bolstered more recently by the discovery of a com- 
plex phase diagram for this system in applied magnetic 
field and pressure. 

§ (Here we take "unconventional" 
to imply the existence of additional broken symmetries 
beyond the usual gauge U(l) broken in classic supercon- 
ductors. 1^) 

Several current Ginzburg-Landau (GL) theories of the 
UPt^ phase diagram attribute the existence of multiple 
superconducting phases to two nearly degenerate super- 
conducting states, either (i) split by a symmetry-breaking 
field, such as the ordered antiferromagnetic moment in 
the basal plane, or (ii) "accidentally" degenerate. |^ Such 
theories can at the same time describe qualitatively the 
anisotropy of the superconducting state, insofar as the 
GL parameters can be chosen to stabilize an order pa- 
rameter at low temperatures and fields which allows for 
a larger number of quasiparticle excitations with wave 
vector in the basal plane. Such a state is strongly indi- 
cated by analyses of both ultrasound jl] and ther- 
mal conductivity measurements. Beyond this crude 
statement, little is known for certain about the exact 
anisotropy or even the symmetry of the superconducting 
state of UPt^. 

Recently, Lussier et al have argued that thermal 
conductivity measurements can shed further light on 
these questions. They showed that the electronic heat 
current dominates the phononic current down to low tem- 
peratures for their high quality samples, and that the re- 
laxation rate I/t^ in the normal state is nearly isotropic. 
Furthermore, their measurements imply the existence of 
large anisotropy in the superconducting state which does 
not simply reflect normal state anisotropy; together with 
transverse ultrasound measurements, these data provide 
convincing evidence for a highly anisotropic gap in UPt^. 

In order to determine the actual gap anistropy for 
U Pt^ , it is necessary to go further and attempt to model 



the data. While it has been stated that such fits to trans- 
port properties cannot be expected to fix the detailed 
anisotropy due to uncertainties in the form of the impu- 
rity scattering amplitude, Fledderjohann and Hirschfeld 
Q argued recently that ratios of transport coefficients 
should lead to more robust conclusions since they can 
depend only weakly on the relaxation times. They there- 
fore focussed on the ratio Kc/nb between the conductiv- 
ities measured for heat currents directed along the c- 
and b-axes, respectively, comparing the data of Lussier 
ct al [Q to weak-coupling BCS calculations using order 
parameters representative of the Eig and i?2« symmetry 
classes of the Dq^ space group of the hexagonal crystal. 
While both states have lines of order parameter nodes 
(and hence higher density of excited quasiparticles) in 
the basal plane, the E2u state has point nodes along the 
c-axis where the order parameter vanishes quadratically, 
in contrast to the linear behavior in the Eig state. In 
consequence, the thermal conductivity (and indeed all 
current-current correlation functions) was found to be 
isotropic in the i?2u state over a spherical Fermi surface, 
despite the intrinsic anisotropy of the superconducting 
state. Ellipsoidal Fermi surfaces do not change the value 
of the normalized conductivity ratio {nd knc) / [^b/ i^Nb) 
from unity in this state ( refers to the normal state) , 
but it is clear that the true hexagonal crystal structure 
will do so. Furthermore, is not clear whether this result 
is specific to the particular i?2u state analyzed, or would 
hold for a more general i?2u state. Understanding the 
extent to which these factors might improve the agree- 
ment with the large measured anisotropy is crucial to the 
E2u scenario proposed by Sauls and Norman |^ in which 
several problems characteristic of GL theories of type (i) 
above are resolved. 

In this paper, we study the influence of both Fermi 
surface anisotropy and gap anisotropy on superconduct- 
ing state transport coefficients, focussing on the thermal 
conductivity data of Lussier et al. [tIJioI] In the first part, 
we use a simple ellipsoidal Fermi surface fit to normal 
state transport data and analyze all gap functions repre- 
sented by ellipsoidal harmonics up through L=5, treating 
various impurity scattering rates, impurity phase shifts. 
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and inelastic scattering effects. First, we find that a fi- 
nite, non-zero T=0 value of the thermal conductivity ra- 
tio, Kc{0)/ Kb{0), of intrinsic origin occurs for a number 
of harmonics, not just those of E2u symmetry. Second, 
we find that the data can be fit reasonably well by gaps 
of both Eig and i?2n symmetry, with the latter fitting 
slightly better than the former, although in neither case 
is a pure harmonic realized. These fits could be differen- 
tiated more clearly by (1) extending the measurements to 
lower temperatures or (2) by increasing or decreasing the 
impurity scattering rate, that is, by analyzing cleaner or 
dirtier samples. Although fits using ellipsoidal harmon- 
ics may be somewhat unrealistic, they allow us to obtain 
some useful analytical results, and determine the quali- 
tative features of order parameter anisotropy with some 
confidence. 

In the second part, we turn to the more general 
case, using the multi-sheeted Fermi surface predicted 
from local density approximation (LDA) calculations ||Tl| ] 
which is in reasonable agreement with deHaas-vanAlphen 
(dHvA) experiments, up to an overall mass renor- 
malization. Two types of gap functions are analyzed: 
Fermi surface harmonic and tight binding (ijjl^ . In 
neither case is an adequate fit found to the data for ei- 
ther E2u or Eig with single basis functions, although one 
of the tight binding gap functions of Eig symmetry has 
some promise. In the Fermi surface harmonic case, this 
poor agreement is due to the large number of nodes these 
functions possess which is unlikely to arise out of any mi- 
croscopic gap equation. In the tight binding case, this is 
likely due to the use of a single basis function. Use of a 
mixed basis set in the tight binding case leads to a good 
correspondence to the data in the Eig case. So far, we 
have not found a comparably good fit for the E2u case. 



I. ORDER PARAMETERS AND FERMI 
SURFACES 

Although a variety of models have been proposed for 
the order parameter of U Pt^, we concentrate here on the 
most popular model, that of a two-dimensional group 
representation. The two variants most commonly ex- 
plored have been the Eig model ||T^ and the E2u model 
For a spherical Fermi surface, the gap function can 
be represented by spherical harmonics. A function of 
Eig symmetry first occurs in the L=2, M=l representa- 
tion (d-wave). The E2u case is more subtle since it is 
an odd parity gap and therefore a pseudo-spin triplet. 

The proposed E2u model assumes that the gap is 
a pure spin triplet with only one component [Sz — 0) 
condensed, however. In this case, £'211 first occurs for 
I32 (f-wave). The £'211 model based on Y32 was originally 
proposed since (i) its nodal structure was similar to 
the previously considered Eig model based on I21, with 
line nodes perpendicular to the c-axis and point nodes 
along the c-axis as indicated by transverse ultrasound 
as well as point contact spectroscopy, ||l^ and (ii) 



it has an upper critical field anisotropy consistent with 
experimental data |po| given the Sz = orientation of 
the triplet order parameter, as demonstrated earlier by 
Choi and Sauls |2l| ] (singlet order parameters give an in- 
correct anisotropy). Sauls |^] in turn showed that this 
model solved a major problem of the previously consid- 
ered Eig model, in that it could explain the existence of a 
tetracritical point in the H-T phase diagram for all orien- 
tations of the magnetic field as observed experimentally, 
at least for axial symmetry. Recently, Park and Joynt l22] 
have proposed that Eig can avoid the problem of an in- 
correct upper critical field anisotropy if the normal state 
Pauli susceptibility has opposite anisotropy to the ob- 
served normal state susceptibility (the latter likely being 
van Vleck dominated). It can also give a phase diagram 
which has a near tetracritical point for certain choices 
of the GL coefficients, with the additional claim that it 
gives a better explanation of the pressure-temperature 
phase diagram than £211 • 

A potential method of resolving these controversies 
would be to obtain more knowledge of the actual form of 
the gap anisotropy. The Eig model has a linear disper- 
sion of the quasiparticle energies about the point nodes, 
whereas the E2U model has a quadratic dispersion. This 
can have a significant effect on transport quantities, as 
pointed out by Yin and Maki. Fleddcrjohann and 
Hirschfeld ||^ exploited this to show that the thermal 
conductivity anisotropy ratio, Kc(0)/Kb(0), is small for 
the Eig case (it would be zero in the clean limit), but is 
unity for the i?2« case, at least for an ellipsoidal Fermi 
surface, with the data of Lussier et al |^ lying between 
these two results but being more consistent with Eig than 
E2u- This in turn motivated Lussier et al to take data 
at lower temperatures, where they conclude that the ex- 
trapolated T=0 anisotropy ratio of about 0.5 is probably 
intrinsic and thus consistent with an E2U model. JTof 

The above analysis of Fleddcrjohann and Hirschfeld j|] 
was based on a particular spherical harmonic form of the 
order parameter on an ellipsoidal Fermi surface. For a 
real metal like UPt^, we would expect that the actual 
order parameter is more complicated, just as we know 
that the actual Fermi surface is multi-sheeted and shows 



strong deviations from axial symmetry. |12| The latter is 
particularly important since the i?2« result Kc/^b = 1 is 
a consequence of axial symmetry. 

To analyze this in more detail, we first consider the 
simple ellipsoidal case treated previously, but look at 
other harmonics besides I21 and ¥32. The conversion 
from spherical to ellipsoidal harmonics can be achieved 
by replacing sin(0) by sin{9) / y'nir and (previously 
unity) by cos^{9) + sm'^{d)/mr HJ] where the mass ra- 
tio, m,. = m±/mc, is equal to 2.8 based on normal state 
transport data. 0| Note that this conversion simply mul- 
tiplies Y21 and 132 by an overall constant, so the re- 
sults for these two cases are independent of the mass 
ratio. This is not true in general. For the Eig case, the 
next higher harmonic to consider is 1^41 ; for the E2U case 
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{Sz = 0), ¥52. Although such higher harmonics seem ex- 
otic, they do play a significant role in certain microscopic 
theories. 

For the real Fermi surface case, we utilize the surface 
obtained from an LDA calculation. This surface, 

which is shown in symmetry planes of the zone in Fig. 1, 
is in reasonable agreement with dHvA data, except for 
mass renormalization effects. The mass renormalization 
would play no role here unless it was anisotropic. Unfor- 
tunately, there is not enough data available to model this 
anisotropy, although the lack of observation of dHvA or- 
bits for fields along the c axis suggests that the renormal- 
ization is anisotropic and will act to increase the mass ra- 
tio, nir. ||l2| To check this, we calculated < v1 > / < > 
where Vb and Vc are the Fermi velocities and <> is an av- 
erage over the Fermi surface determined by using a linear 
tetrahedron decomposition of the Brillouin zone. [ ^J28| 
This quantity, equal to the ratio uncI i^-Nb given the ob- 
servation of isotropic relaxation time, Q is 2.8 from ther- 
mal conductivity and 2.7 from resistivity. The LDA 
calculation gives 2.1. This is consistent with the above 
observation that the mass renormalization anisotropy 
acts to increase m^. Since we are comparing in this paper 
to data normalized to the normal state value, we ignore 
this mass ratio discrepancy in this paper, although we 
caution that it could influence some of the results pre- 
sented here. In particular, an alteration in Fermi surface 
topology would certainly change the gap anisotropy, but 
even a momentum dependence of the mass renormaliza- 
tion would affect the results since this would act to alter 
the weighting of various momentum vectors in the equa- 
tions. 

The problem of what gap function to use for the real 
Fermi surface case is a more complicated issue. Allen 
showed a number of years ago that functions could be 
constructed on the Fermi surface which were orthornor- 
mal and so could be used as basis functions, which he 
labeled Fermi surface harmonics. [T^ ] They can be ob- 
tained from the spherical harmonics by replacing ki by 
Vi, [ p4[ although there are more Fermi surface harmon- 
ics for a given power of Vi than there are spherical har- 
monics since the lattice symmetry is discrete rather than 
continuous. A complication is that the relative size of 
the gap function on independent sheets of the Fermi sur- 
face cannot be determined outside the context of a mi- 
croscopic theory of the superconductivity. In fact, the 
problem is a more general one since the simple gap func- 
tions treated here are likely to be modulated at each k 
vector by some complicated function (of Aig symmetry) 
that contains the affects of wavefunction anisotropy, etc. 
We further note that the Fermi surface harmonics will be 
small when the Fermi velocities are small. On the other 
hand, we would anticipate that in most microscopic the- 
ories, the gap function will be large where the f electron 
weight is highest, which is also where the Fermi velocities 
will in general be lowest. Keeping these caveats in mind, 
we define Fermi surface harmonics for the Eig case and 
case by replacing ki by Vi in the I21 and I32 spher- 
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VcVr and the E2U gap VcV^ where — + v^. |^ The 
nodal structure of these functions are very complicated 
given the complicated Fermi surface geometry. There are 
many points on the Fermi surface where the velocity vec- 
tor points either along or perpendicular to the c-axis. In 
any of these cases, the Eig and i?2u Fermi surface har- 
monics will vanish. 

An alternate set of basis functions can be generated 
by tight binding expansion. In the square lattice case, 
the lattice vectors of type (1,0) lead to a d-wave state of 
the form cos(fca;) — cos(fcj,), which is currently the lead- 
ing model being explored for high temperature cuprates. 
For the hexagonal closed packed case, Eig and E-2u first 
appear for primitive lattice vectors of the type (0,1,1). 
These can be generated from the next near neighbor ba- 
sis functions listed by Putikka and Joynt 10 by multi- 
plying their Eiu and E2g functions by sin(fc2c). They are 
for Eig 

1 %/3 

/i = \/2sin(fc2c) cos(-fcya) sin(— fcj^a) 

2 1 a/3 

/2 = sin(fc2c)(sin(fcj^a) -f su\{-kya) cos{—k^a)) (1) 

and for E2U 

fi = -^sin{kzc){cos{kya) - cos{^ky a) cos{^ k^a)) 

1 \/3 
/2 = V2sm{kzc) sm{-kya) sin( — fcj;a) (2) 



with a gap modulus of y/ fi + fi ioi an assumed gap of 
the form /i +if2 (the l,i state). Both functions have line 
nodes in the k^ — and k^ = tt/c planes. In addition, 
the Eig function has point nodes with linear dispersion 
along all three symmetry axes {T — A, AI — L, K — H), 
whereas E2U has quadratic point nodes along F — ^ and 
linear point nodes along K — H . We note that the ba- 
sis functions listed by Putikka and Joynt for the near 
neighbor case are not properly invariant under recipro- 
cal lattice translations due to the non-symmorphic na- 
ture of the UPt^ lattice (that is, the near neighbors are 
separated by a non-primitive translation vector). This 
problem has been addressed by Konno and Ueda. |l^] 
Under the highly simplistic assumption that the phase of 
the single particle wavefunctions on the two sites in the 
unit cell is determined by a simple near neighbor inter- 
action, they were able to generate analytic near neighbor 
basis functions which have proper translational symme- 
try. In this paper, we use their F^ {Eig) and {E2u) 
basis functions. These, in fact, can be generated from 
the tight binding basis functions discussed above by re- 
placing sin(fczc) by sm(kzc/2) where (j) is the Fourier 
transform of the near neighbor distance vectors projected 
onto the basal plane 
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((/) is complex since the lattice is non-symmorphic). The 
effect of this is to remove the line nodes in the kz — tt/c 
plane and the linear point nodes along K — H. 



II. THERMAL CONDUCTIVITY 

The thermal conductivity k in the presence of impu- 
rities is evaluated using a Kubo formula for the heat- 
current correlation function as in the original treat- 
ment for an s-wave superconductor by Ambegaokar and 
Tewordt. ^ This treatment was generalized to uncon- 
ventional states by several groups, giving results for a 
spherical Fermi surface and model p- and d-wave states 
which agreed qualitatively with experiment. fl,|5|, ^ , ^ 
In the limit of vanishing impurity concentrations, iden- 
tical results were also obtained by Arfi et al Q using a 
transport equation method. 

The Kubo formula approach begins with an impurity- 
averaged single-particle matrix propagator 
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where r* represent the Pauli matrices spanning particle- 
hole space. Here, we have already exploited the ap- 
proximate particle-hole symmetry of the normal state, 
as well as the symmetries of the gap functions which 
lead to vanishing off-diagonal scattering self-energy con- 
tributions. In this limit, only self-energy contibutions 
to the frequency cj, namely a) = cj — Sq need to be in- 
cluded. [gH The self-energy Sq due to the elastic im- 
purity scattering is treated in a self-consistent T-matrix 
approximation and is given by So = rGo/(c^ — Gq), 
where F = nin/ {'kN^) is the unitarity limit scatter- 
ing rate depending on the concentration of defects n^, 
the electron density n, and the density of states at the 
Fermi level Nq . The quantity c = cot 5q parameter- 
izes the scattering strength of an individual impurity 
through the s-wave phase shift Sq. In this work we con- 
sider primarily unitarity limit scattering c = since it 
is clear that weak scattering will lead to a weak temper- 
ature dependence inconsistent with experiment for the 
states in question. [p4|j35| The integrated propagator is 



Go = (l/27riVo)EkTr{T^5(k,tj)}. The equation for the 
self-energies are now solved self-consistently. We ignore 
the complication of resolving the gap equation in this 
paper since all results are scaled to Tc and impurity cor- 
rections are small. The gap itself is given by 
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where A bcs (T) is the BCS gap and /^r is one of the basis 
functions (or mixtures thereof) described in the previous 
section. 

The bare heat current response is given by a convolu- 
tion of the Green's function g with itself at zero exter- 
nal frequency and wave vector weighted with the bare 



heat current vertex luvi^t^. Impurity scattering ver- 
tex corrections to current-current correlation functions 
have been shown to vanish identically for even parity 
states (Ak = A_k)- Even for odd parity states, such 
corrections vanish in the unitarity limit. For the diagonal 
thermal conductivity tensor one obtains 
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where cD' and lu" are the real and imaginary parts of uj 
and Fjv = F/(l-t-c^)is the normal state scattering rate. 
Here dSk is the area measure on the Fermi surface, and 
Wk is the Fermi velocity. 

For a complete description of the data, we must take 
into account the effects of inelastic scattering. This is 
known to vary as bT^ times the elastic rate in the nor- 
mal state, with b ~ A/K'^. This effect can be included 
in the above equations by replacing F by r(l -f bT^). 
In the superconducting state, we can make the ansatz 
that the inelastic rate varies as bT^/Tc since the num- 
ber of quasiparticles varies as T/Tc at low temperatures 
due to the line nodes in the gap. The exact form of 
this makes little difference, since by far the largest effect 
inelastic scattering has is on kn{T). Therefore, for prac- 
tical purposes when comparing to data normalized to its 
value at Tc, one can simply scale the result of Eq. 6 by 
(l + 6T<?)/(l-h6T2). 

Finally, we note that UPt^ has a split superconducting 
phase transition. In the E models considered here, this 
splitting is assumed to be due to the weak antifcrromag- 
netism which has orthorhombic symmetry. Its effect is 
to cause only one of the two E components to condense 
at the upper phase transition. Thus, in the region be- 
tween the upper {Tc+ = 0.50K) and lower (rc_=0.44 K) 
transitions, the point nodes along the c axis become line 
nodes perpendicular to the basal plane. This explains 
the lack of gap anisotropy in the thermal conductivity 
observed in this region. Below Tc-, the second E com- 
ponent condenses, and the anisotropy begins to occur, as 
observed. Modeling this is complicated since the calcula- 
tion would have to be performed for orthorhombic sym- 
metry with two different gaps and appropriate domain 
averaging performed. We therefore take the approach of 
previous work which ignores this symmetry breaking but 
normalizes Kc/ Hb to its value at Tc- • This normal- 

ization does not work so well, though, when comparing 
to the individual Ki{T)/T themselves, since the thermal 
conductivity does change above Tc- ■ We have found that 
normalizing Ki{T)/T to its value at Tco = 0A7K, the av- 
erage of Tc+ and Tc- (the "hexagonal" Tc), works quite 
well in this regard. Obviously, one cannot take too se- 
riously the results in the immediate vicinity of Tc until 
the effects of the symmetry breaking field are properly 
included. 
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III. RESULTS 

A. Ellipsoidal harmonics 

Results for Kc/ Hb are shown in Fig. 2 for all harmonics 
through L=5 on an ellipsoidal Fermi surface with = 
2.8, compared to the experimental results of Lussier et 
al. [0 This quantity is normalized to its value at Te- 
as discussed above. The results were generated with an 
impurity scattering rate in the unitarity limit of 0.1 Tc, 
consistent with experimental data (particularly with the 
observation of a residual linear specific heat coefficient of 
0.16 the normal state value). |]l^ To understand these 
results more clearly, we have analytically calculated 

in the clean limit for a spherical Fermi surface, and show 
these results in Table 1. We see that harmonics of the 
form 111, which have only point nodes along the c axis, 
give a divergent ratio. On the other hand, only two of 
the remaining harmonics, Yio and ^21, give a ratio of zero 
(the non-zero value in Fig. 2a is due to impurity-induced 
gaplessness ) , with the rest giving an intrinsic non-zero 
ratio. In particular, we note that harmonics of the form 
Ilo only have line nodes, so a finite, non-zero ratio is not 
something just associated with quadratic point nodes or 
with gaps of E2u symmetry. We also note that no pure 
harmonic provides a good fit to the observed anisotropy, 
Q although some higher harmonics give adequate fits. 

To study this further, we have looked into the possi- 
bility of mixed solutions. For Eig, we included mixing 
of 1^21 with 141; for E2u, F32 with ¥52- Typical re- 
sults for Kc/ Kb are shown in Fig. 3a, with the coefficients 
roughly optimized to fit the data. Both fits give a reason- 
able description of the data, with the lowest temperature 
data intermediate between the two results (it should be 
remarked that the error bars on the experimental nd Hb 
are about 15% at low temperature [|l^). The Eig fit can 
be greatly improved at lower temperature by going to a 
larger scattering rate of 0.3 Tc, but the individual Ki/T in 
this case are in poor agreement with experiment. Alter- 
ing the scattering phase shift from the unitarity value of 
cot((5o) of to 0.2 slightly improves things at the lowest 
temperatures, but this is probably not significant given 
the experimental error bars. For the E2U case, lowering 
the scattering rate by a factor of ten only slightly sup- 
presses the ratio and only for temperatures below where 
experimental data exist (the same slight suppression also 
occurs by increasing cot((5o)). In Fig. 3b, we compare 
these fits to the individual Ki/T, normalized to their 
value at Tco as discussed previously. As can be seen, 
both Eig and i?2u provide good fits to the data, with 
the E2U fit being slightly superior. An interesting point 
is that the experimental Ki/T are linear in temperature 
down to the lowest measured temperatures. This behav- 
ior cannot continue indefinitely since ki/T would be zero 



at a temperature larger than zero (that is, ki/T must 
fiatten off). The calculated curves, though, predict that 
this flattening should occur in the measured temperature 
range in contradiction with experiment. As suggested by 
the authors of Rcf . jl0| , the calculated low temperature 
behavior can be improved by reducing F to roughly one 
tenth its normal state value. Although this does improve 
the fit at low temperatures, it leads to a substantial devi- 
ation from the data at higher temperatures. In fact, we 
have found that the value V /Tc ~ 0.1 (the normal state 
value) gives roughly the best fit over the entire temper- 
ature range below Tc- We have also found that altering 
the scattering phase shift from the unitarity value of 7r/2 
does not improve the fit in this regard, at least for small 
values of cot((5o). If this discrepancy between the low 
temperature and high temperature behavior is taken at 
face value, a strong temperature dependence of either 
c = cot((5o) or r must be assumed. Although dynamical 
scattering effects could easily influence the phase shifts 
in this way, we have been unable to flnd a satisfactory 
phenomenological explanation of the data in these terms. 
While a T-dependent parametrization of F could possibly 
account for the discrepancy, we have no physical under- 
standing of how such effects could arise. 

The gap functions in the Eig and E2U cases are plotted 
as a function of polar angle in Fig. 4. Both gap functions 
look similar (except for the different dispersions around 
the point nodes at zero degrees). This indicates that the 
primary determinant of the thermal conductivity is the 
overall shape of the gap function. We also note that the 
maximum gap occurs at a polar angle of 52 degrees for 
Eig and 49 degrees for E2u- Not only are these angles 
close, but they are also close to the angle of 54 degrees 
that the vector connecting near neighbor uranium atoms 
makes with the c axis. This could be taken as indirect 
evidence that the electrons in the Cooper pairs reside at 
near neighbor sites as would be predicted by microscopic 
models based on antiferromagnetic spin fluctuations. 

We conclude this part by remarking that both the Eig 
and E2U models can explain the data. One way to more 
clearly distinguish between the two would be to carry 
the experiments to lower temperatures, although given 
experimental error bars, it may be difficult to conclude 
anything definitive. At the least, one would hope to see 
Ki/T flatten off. Perhaps a better way would be to de- 
grade the quality of the sample. In the presence of a finite 
concentration of impurities, states with line nodes yield 
a linear term in the thermal conductivity, Ki{T) ^ T/Aq 
at the lowest temperatures. For generic configurations of 
the thermal current jg and the line nodes, the propor- 
tionality constant is actually independent of the impurity 
scattering rate to leading order, yielding a universal value 
for the low-T thermal conductivity analogous to the elec- 
trical conductivity result found by P.A. Lee. In the 
final stages of writing, we received a paper from Graf et 
al. [ |39|] in which this result was obtained independently 
and explored in some detail. 
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Because of impurity-induced gapless effects of this 
type, the anisotropy ratio Kc(0)/k6(0) is always finite 
even for states hke -Eig, as mentioned above. In such 
a situation, an estimate of the anisotropy may be per- 
formed by considering Eq. 7 in the gapless regime, i.e. 
take uj — auj + ij, with a, 7 constant. We then find 
K40)/n,{0) = {v^,Fi)/{vlFi), with F, = (7^ + Al)-^/\ 
For the Y21 case (Eig) in spherical symmetry, we find 
Kb{T) ~ T/Ao, ~ 7T/AI giving k,(0)/k6(0) ~ 

27/ Aq for small concentrations. The residual broaden- 
ing 7 is found by solving the transcendental equation 

+ 7^(-Fa:)^ = r(i^fc), and yields a square root depen- 
dence on concentration, Kc{0) / Kb{0) ~ (F/Ao)^/^ up to 
logarithmic corrections in the unitarity limit c = 0. In 
Fig. 5, we plot the impurity concentration dependence 
of the anisotropy ratio in both Eig and cases for 
T=0.02 and 0.2 Tc- We note the qualitatively stronger 
dependence of the anisotropy ratio on impurity scatter- 
ing rate for the Eig case as compared to the E2U case at 
low temperatures. 

B. Fermi surface harmonics 

The above analysis assumes an ellipsoidal Fermi sur- 
face. The actual Fermi surface for UPts is very compli- 
cated and could substantially alter the above conclusions. 
In this part, we present results using Fermi surface har- 
monic gap functions as described in the first section with 
the Fermi surface shown in Fig. 1. We again assume 
an impurity scattering rate of 0.1 Tc- Results for Kc/Kt 
are shown in Fig. 6a and for Ki/T in Fig. 6b. The Eig 
case gives a fair representation of the anisotropy ratio, 
especially at higher temperatures. Both cases, though, 
predict too large Ki/T at lower temperatures. This oc- 
curs since the Fermi surface harmonic gap functions have 
a very large number of nodes. These nodes are due to 
the complicated Fermi surface, which has many points 
on it where the velocity vector is either parallel or per- 
pendicular to the c axis, either case in which the Eig and 
E2U Fermi surface harmonics vanish. This large number 
of nodes is unlikely to arise out of any microscopic gap 
equation since such a solution would not have an optimal 
condensation energy. This indicates that Fermi surface 
harmonics are unlikely to be useful in modeling heavy 
fermion superconductors. Because of this, and since the 
addition of higher order harmonics will not reduce the 
number of nodes, we have not explored using mixed Fermi 
surface harmonics as done in the previous subsection for 
ellipsoidal harmonics. 

C. Tight binding functions 

We next present results for tight binding gap func- 
tions. These functions represent short range interactions 
in the lattice, and therefore are more likely to arise out of 



a microscopic gap equation than the Fermi surface har- 
monics. Results for the functions based on (0,1,1) type 
primitive lattice vectors are shown for Kc/ Hb in Fig. 7a 
and Ki/T in Fig. 7b. Although the magnitude of Ki/T 
is improved over the Fermi surface harmonic case, the 
anisotropy ratio is in poor agreement with experiment. 
We therefore turn to results using the Konno-Ueda func- 
tions based on near neighbor interactions |l^ (Figs. 8a 
and 8b). The observed anisotropy ratio is intermediate 
between the Eig and E2U cases. Although Ki/T is too 
large at lower temperatures for E2m it is not too bad for 
Eig. Perfect agreement with experiment would not be 
expected anyway since the above basis functions do not 
take into account the complicated single particle wave- 
functions which occur in heavy fermions due to f orbital 
degeneracy and f-ligand hybridization. It is interesting 
to note that the near neighbor tight binding functions 
provide the best overall comparison to the data for those 
functions we have analyzed on the real Fermi surface, 
since this reinforces the idea of near neighbor pairing that 
was suggested from the ellipsoidal results above. Given 
this, it will interesting in the future to calculate k for re- 
cent microscopic models which do not involve such pair- 
ing. g 

As in the ellipsoidal harmonic case, to test the idea of 
whether a mixture of tight binding basis functions would 
improve the results, we have done calculations mixing 
the two tight binding functions considered above for each 
symmetry. Only a rough optimization of the mixing co- 
efficients could be determined due to calculational de- 
mands. These results are presented in Figs. 9a and 9b. 
A good fit was obtained for the Eig case with an equal 
admixture of the two functions. This result is some- 
what surprising since an equal admixture would imply 
(within a tight binding framework) that the pair inter- 
action has a substantial range (the near neighbor sepa- 
ration is 7.8 a.u., but an (0,1,1) vector is 14.3 a.u.). In 
the E2U case, we never found an adequate fit, although 
we do show a typical result (equal admixture, but with 
opposite sign). One should be cautious, though, about 
ruling against an E2U model based on this, since there are 
other functions involving non-primitive translation vec- 
tors, with distances comparable to the (0,1,1) primitive 
vectors, which we have not considered here (one at 13.4 
a.u., another at 15.2 a.u.) 

IV. CONCLUSIONS 

In conclusion, we have found by analyzing simple ellip- 
soidal models for the Fermi surface that recent thermal 
conductivity data cannot unambiguously differentiate be- 
tween the Eiu and E2U models for the symmetry of the 
order parameter. Such a differentiation should be possi- 
ble by looking at samples with degraded quality, since in 
the E2U case, the non-zero value of Kc(0)/Kb(0) is intrin- 
sic, whereas for Eig, it is due to impurities. We have also 
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found that results using realistic Fermi surfaces and gaps 
with proper lattice translational symmetry differ signifi- 
cantly from those based on ellipsoidal Fermi surfaces, and 
have discovered an Eig tight binding gap function which 
gives a good representation of the data. We emphasize 
that the thermal conductivity data appear to put great 
constraints on the overall shape of the gap function, and 
thus will be an important ingredient in determining the 
validity of microscopic models for the superconductivity 
in UPta. 
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TABLE I. Kc(0)/k(,(0) in the clean limit for several Ylm 
gap functions for two values of the mass ratio, m^. The har- 
monics Y21 and 1^41 have Eig symmetry and the harmonics 
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FIG. 1. LDA Fermi surface for UPt^ plotted in the symme- 
try planes of the Brillouin zone. The surface is composed of 
five bands, three centered around F and two centered around 
A. 

FIG. 2. hic/i^b (normalized to its value at Tc-) for various 
harmonics through L=5 (curves labled by L,M) on an ellip- 
soid (mr=2.8). All plots, unless otherwise noted, are for an 
impurity scattering rate, F, of 0.1 Tc in the unitarity limit, 
with an inelastic scattering rate 4r^ times the elastic rate. 
The black dots are data from Ref. Urn. 



FIG. 3. (a) hie/ lily and (b) k^/T for mixed harmonics on 
an ellipsoid (mr=2.8). The curve marked Ig is for E\g 
{Y21 - 0.15^41) and the one marked 2u for E2n^ ("^32 +0.2y52). 
Black dots in (a) and black (white) dots in (b) are data from 
Ref. 0. The solid (dashed) curves in (b) are the respective 
theoretical results for i=b (i=c). 

FIG. 4. Order parameter versus polar angle for mixed har- 
monics on an ellipsoid (mr=2.8). 



FIG. 5. Kc/k6 at r/Tc = 0.02 (solid line) and T/T^ = 0.2 
(dashed line) vs. normalized impurity scattering rate, F/Tc, 
for pure ellipsoidal harmonics ¥32 (2u) and y2i (Ig), as well 
as for the mixed harmonics of Fig. 3a (2u* and Ig*). 



FIG. 6. (a) Kc/k6 and (b) KijT for Fermi surface harmonics 
on the LDA Fermi surface. Same notation as Fig. 3. 

FIG. 7. (a) nc/nb and (b) K,i/T for tight binding functions 
(0,1,1 type primitive lattice vectors) on the LDA Fermi sur- 
face. Same notation as Fig. 3. 

FIG. 8. (a) Kc/Kb and (b) Ki/T for near neighbor tight 
binding functions (Konno-Ueda type) on the LDA Fermi sur- 
face. Same notation as Fig. 3. 

FIG. 9. (a) Kc/Kb and (b) Ki/T for mixed tight binding 
functions on the LDA Fermi surface. Same notation as Fig. 
3. 
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